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Abstract 

In this paper, two classes of exact analytic time-dependent soultion of magnetic annihi- 
lation for incompressible magnetic fluid, have been obtained by solving the magnetohydro- 
dynamic (MHD) equations directly. The solutions derived here possess scaling property with 
time t as the scale factor. Based on these two solutions, we find that, for some given inflow 
fields, the evolution of the annihilating magnetic field can be described by the solutions of 
certain ordinary differential equations whose variables are dilated simply by time t. The 
relevant evolution characteristics in the process of magnetic annihilation are also revealed. 

PACS number(s): 95.30.Q, 52.30, 75.50.M 



As is well known, magnetic reconnection is a very important physical phenomena in 
many astrophysical objects. It is believed that, magnetic reconnection can serve as such a 
mechanism to transfer the magnetic energy into heat and kinetic energy of the plasma Q . 

Generally, due to the nonlinearity of the MHD equations, it is very difficult to find 
its analytic solutions. We have to turn to the numerical simulations. But, unfortunately, 
numerical simulations are usually limited by artificial resistivity and the implementation of 
sensible boundary conditions ||. [| . 

Anyway, the exact analytic solutions are needed and helpful. On one hand, they can 
provide the important insights into the problem concerned, and on the other hand, they can 
be used to test the relevant numerical methods. Just based on the beautiful exact solutions, 
Parker found that, diffusion is an essential part of the reconnection of lines of force at a neutral 
point ||. Since then, a lot of excellent analytic solutions have been found and applied to 
analyze some concrete models [j], ||, ||, 10]. 



Recently, Watson, Craig and Henton etc. have made a great progress in searching for the 
exact solutions []TT| , 12, 13, They found that, exact families of reconnection solutions can 



be constructed in both two and three dimensions, and, these solutions display the essential 
characteristics of fast reconnection. Furthermore, they have shown that, the reconnection 
solutions can be obtained by superposing the velocity and magnetic fields of simple magnetic 
annihilation model. 

So, annihilation model can provide a solid basis for our understanding the reconnection 



process. Up to now, there are few works on the time-dependent case [|H], 15 1. Most of the 
previous works in searching for the exact solutions concerntrate on the case of the steady 
MHD. But, it is desired to remove the restrictions of stationarity and incompressiblety etc. 
for the realistic application. 

In this paper, we present two time-dependent analytic solutions of the incompressible 
MHD. We begin by introducing the basic equations in section 1. The solutions are obtained 
in section 2. And, the relevant physics indicated by the solutions is analysized in section 3. 
Finally, the conclusions are summarized in section 4. 

1. The Basic Equations 

The incompressible MHD equations can be written as the following nondimensional form 



0: 



dv 



-(v ■ V)v- VP + (V x B) x B (1) 



3B 

— = V x (v x B) + i]V 2 B (2) 

V-B = (3) 

V-v = (4) 
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where the length, magnetic field B, velocity v, time t, gas pressure P and the resistivity rj have 
been non-dimensionalized by the typical coronal parameters - length-scale L c , background 
field strength B c , Alfven speed va = B c /(8irp) 1 ^ 2 , Alfven travel time ta = L c /va, magnetic 
pressure B 2 /(8ir) and vaL c respectively. 

By taking the curl of Eq.(|l|), we can remove P in above equations and get 

^^ = -Vx[(fV)»1+Vx[(Vx5)xB] (5) 

Now, we need to deal with Eqs.(|2|, [3l ||, ||). Once we solve them and get v and B, we can 
calculate the pressure P according to Eq. (]!]). 

2. Magnetic Annihilation Solution 

In order to solve Eqs.(|2[ || |||, |B|), suppose that B, v have the following forms 

B = {0, Y(x,t), 0} (6) 
v = {U(x,t), yV(x,t), 0} (7) 

for the annihilation solution. Here, U,V,Y are the functions of x and t. And thus, Eq.(^) is 
satisfied automatically. 

Now, Eqs.(|, |) and Eq.(|) give 

—U(x,t)+V(x,t) = (8) 

—V{x,t) + U(x,t)— V(x,t) + V(x,t) 2 + f(t) = (9) 

-Y(x,t) + —(U(x,t)Y(x,t))-r l ^Y(x,t) = (10) 

where, f(t) is a function of t to be determined later by the boundary condition. 
If we further suppose that 

U(x,t) = U(x)g(t) 

from Eq.(P), it is evident that 

V(x,t) = V(x)g(t). 
Substitute the above two equations into Eqs.(||, |9|), we achieve 

^-g(t) + C ig (t) 2 = (11) 

f(t) + c 2 g(t) 2 = (12) 
^U{x) + V{x) = (13) 

U{x)^-V{x) + V{xf - Cl V(x) -c 2 = (14) 
ox 
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with Cj (i = 1,2,3, ...) being the (integral) constants. Now, we are in position to state how 
to solve g(t), U(x) and V(x) from the above equations. 
Eq.fljl]) tells us that 

9(t) 



1 



Starting out from Eq.(13, 14), we can obtain two solutions. 
A. Take U(x) as the linear function of x, then we have 



U(x) 
V(x) 

C2 



C 4 X + c 5 
-c 4 

C1C4 + c\. 



B. Let U(x) = P(x) + C4, and substitute it into Eqs.([D| 14), we arrive at 
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dx 



P{x) + V(x) = 



P(x)— ^(x) + V(x) 2 + C4^^(x) - ciV(x) - c 2 



0. 



(15) 



(16) 
(17) 
(18) 



(19) 
(20) 



if we further have P(x) = c^V(x) and ^V(x) = C6^(x), Eqs.(19, ^) can be satisfied easily. 
Then, we get 



V(x) 
U{x) 

CQC4 
C5C6 
C-2 



c 7 e C(iX 

c 5 c 7 e C6X + c 4 

ci 

-1 

0. 



(21) 
(22) 
(23) 
(24) 
(25) 



In the next step, inserting these solutions into Eq.(p^l), we can finally solve the magnetic 



field. 

To give a feeling on the form of the magnetic field B, let us consider the asymptotic 
solution of Y(x, t). When 77 — > 0, Eq.(|l0|) is simplified as 



(26) 



-Y(x,t) + —(U(x,t)Y(x,t)) = 0. 



Take Eqs.(15, 16) and Eq.(22) into account, we get 



Y(x,t) = c 9 



(cit + C 3 ) c l 

{C4X + C 5 ) c 4 



for solution A, and 



(27) 



Y(x,t) = cge C6X (c4e c&x + C5C7) c ^ (cit + c 3 ) c i 



(28) 
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for solution B. 

With a suitable choice of the relevant constants, we can obtain 

Yix - ( > - l^fTT? <29> 
nx - t] ~ FWT? <30> 

corresponding to A and B. Here, the singularity exists evidently at the neutral interface x = 
for both A and B. The presence of the singularity is naturally unavoidable when diffusion 
vanishes, as discussed by Parker H. 

Eq.([U]) tells us taht, Y(x,t) can not be simply decomposed as the form of X(x)T(t), 
after the similar supposition for U(x, t) is made. To solve the magnetic field B analytically, 
a certain transformation is needed. 

For solution A, U(x,t) = §Jf^. Suppose Y(x,t) = Y(s), and 

a = ( C4 x + c 5 )( Cl t + c 3 ) a (31) 
with a = -i, Eq.© can be transformed into 

and for solution B, U(x,t) = C5C J.^ C ^ C4 , similarly, let 

s = c 5 c 7 e cex ic 1 t + c 3 ) a (33) 
with a = — 1, Eq.(jl5) can be read as 

Ws-g^Yis) + (ryce - s)—Y(s) - Y(s) = 0. (34) 

Now, the partial differential equation Eq.(|io|) has been transformed into the ordinary 
differential equations, whose solutions can be expressed as special functions. The solutions 
of Eq.(B3) and Eq.(M) are 



w \ n 17 r c 4 1 (-ci + 2c 4 )s 2 A c 4 3 (-ci + 2c 4 )s 2 

Y{S) = CllFl[ ~^-^> 2' 44v ] + 21 " ^2^' 2' 4cl, ] (35) 

and 

y( s ) = C ie ^ + C 2 e^r[0, — ] (36) 

respectively. Here, C\,Ci are the integral constants, 1F1 is the Kummer confluent hyperge- 
ometric function while T the incomplete gamma function. 

Up to now, we have obtained two classes of the exact analytic solutions for both velocity 
field and magnetic field. We would like to give a short comment here. It is deserved to note 
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that, Eqs.([3l], 33) are actually the scale transformation with time t as the scaling factor. So, 
our solutions will perform certain scaling characteristic, i.e., when the velocity fields are taken 



together with 



as Eqs.fllq, |22| , |15| ), the evolution of the magnetic field with time decribed by the solution of 
Eq.([l(]), can be described by the solution of the ordinary differential equation Eqs.(32, 34) 
dilated simply by time t according to (c\t + c^) a . 

3. Physical Characteristics of the Solution 

In this section, we discuss the physics implied by the above solutions. In doing so, we 
have to fix the free parameters. The relevant functions are taken as 

9(t) = ^ (37) 

U(x) = -x (38) 
V{x) = 1 (39) 

U(x) = e~ x -l (40) 
V(x) = e~ x (41) 

for solution B. In fact, there is another possible choice such as U(x) = 1 — e x , V{x) = e x 
with g(t) given by Eq.(^7|) for the flow field of the annihilation solution, but it is very similiar 
to the one given by Eq.(4Cj4"l|), so we will not discuss it here. 



for solution A, and 



The velocity field for solutions A and B are 

v = { — , 0} (42) 

and 

v = { 6 X ~\ ^4, 0} (43) 

respectively. They are symmetric (for A) or asymmetric (for B) stagnation-point flow weak- 
ening with time. For asymmetric case - solution B, the flow field is shown in Fig.l. 

With the given velocity fields U(x,t), the equation of magnetic field B(x,t) can be read 

as 

7? £ y(s) + r^ y(s) + y(s) = ° (44) 

with the solution 

y(.)—c-*41-%] (— ^t) («) 

and 

VS-Q^Y{s) + {s + ri)—Y(s) + Y(s) = (46) 
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with the solution 

Y(s) = Ce-^m ~] ~ T[0, --]) (s = ^r) (47) 

for A and B correspondingly. To determine the relevant constants, l A | a; =o, t=o = and 
f|| 2=0, t=o = C are required here. 

With the known magnetic fields, we can now solve the current and the energy dissipation 



rate by Joule heat. Based on Eqs.(45, [47]), we achieve 



, . ,1 .5 3 3x .. 5x r 1 1 5 3x ^ . . 

jM = ° { Vm^ 2' -4^TI) ] " 6,(t + l)3/ 2 lFl[ ¥' 2' -4^TI) ]} (48) 

for solution A, and 

^ " - c < 1 + ;*rrr *^< r [°- "iri) 1 - r|0 ' - (49) 

for solution B. Finally, according to 

W„ =< m 2 >= V [ j(x, tfdV (50) 



we can calculate the power provided by Ohmic dissipation. 

The dependence of B(x,t),j(x,i) on (x,t) and W n on time t are shown in Figs. (2-7). 

It can be seen from the figures that, for solution A, the maximum of the magnetic field 
departs from the initial magnetic neutral interface with time. So, the dissipation area will 
get gradually wide. Following it, the current becomes smaller and smaller. The biggest rate 
of the diffuse of the magnetic energy due to the Joule heat takes place at the beginning. 
The solution B is interesting, its current is evidently a soliton. Hence, when the wave-packet 
of the current is located in the domain that we consider, the dissipation power keeps as a 
constant. But, once it passes through the domain, the dissipation power will descend with 
time. 

4. Conclusions and Discussions 

In conclusion, we have obtained two classes of full analytic time-dependent solutions for 
magnetic annihilation by directly solving the MHD equations. 

The solutions which we get here have such a scaling property that the evolution of the 
relevant physical observables are described by the solutions of certain ordinary differential 
equations, with the variables being some functions of x dilated by time t. Besides, in the case 
of asymmetric inflow, a soliton-like current can be formed. Finally, the solutions are used to 
reveal the relevant physical characteristics in the process of the magnetic annihilation. 

However, the time-dependent solutions obtained here are the relative simple ones. Based 
on the method provided by Watson, Craig and Henton etc. [11, 12, 14], we expect to 
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construct the time-dependent reconnection solutions. The work along this direction is under 
way. 
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Figure 1: Stream line for the asymmetric flow field. Where, v(x,t) = j 6 ^ 1 , ^pp, 0}. 
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Figure 2: Magnetic field Y(x,t) versus x in different time t. Where, rj = 0.01, U(x,t) = 

X 

t+1- 
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Figure 3: Current j{x,t) versus x in different time t. Where, rj = 0.001, U(x,t) = — jrrr- 
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Figure 4: Dissipation power of the magnetic energy versus time t in the calculated domain. 
Where, r? = 0.01, -1 < x < 1, U(x,t) = -^. 
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Figure 5: Magnetic field Y(x, t) versus x in different time t. Where, r] = 0.1, U(x, t) 



t+i 
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Figure 6: Current j(x,t) versus x in different time t. Where, ij = 0.1, U(x,t) = e t+1 ' 
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Figure 7: Dissipation power of the magnetic energy versus time t in the calculated domain 
versus time t. Where, 77 = 0.1, — 5 < x < 5, U(x,t) = e t ^~ 1 . 
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